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ABSTRACT 

The basic parameters determining the flow pattern for a nonviscous, non- heat- 
conducting gas in a semidetached binary system without a magnetic field are identified. 
Three-dimensional gas-dynamical modeling of the mass transfer enables investigation 
of the influence of these parameters on the structure of gas flows. The parameter on 
which the flow pattern depends most strongly is the adiabatic index 7. The effect of 
other parameters is small, and only leads to unimportant quantitative changes in the so- 
lutions obtained. The main properties of flows typical of semidetached binaries without 
magnetic fields are summarized. 
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1. Introduction 

Semidetached binary systems belong to the class of 
interacting stars, in which one of the components fills 
its Roche lobe and there is mass exchange between 
the components through the inner Lagrange point. 
Beginning with the pioneering work of Prendergast 
[1] , the gas dynamics of matter flows in semidetached 
systems have been extensively investigated by many 
authors. As a rule, two-dimensional (2D) numerical 
models have been used in calculations (see, for exam- 
ple, [2-9]). Unfortunately, the limitations of the 2D 
approach cast doubt on the adequacy of some results 
obtained using these models. Analytical studies of 
the applicability of 2D models for accretion disks [10] 
have shown that these models give correct results only 
in certain limited cases (with an adiabatic index 7 = 1 
for disks in an external gravitational field and 7 = 2 
for self-gravitating disks). In all other three- 
dimensional (3D) treatment must be employed. The 
comparison [11] of numerical calculations performed 
using 2D and 3D models confirmed the limited appli- 
cability of 2D models, and demonstrated the necessity 
of using 3D models in numerical analyses of the flow 
structures in semidetached binaries. 

Three-dimensional modeling calls for substantial 
computational resources, and for this reason, such 
studies are few in number (the best known are [12- 
21]). In addition, the morphology of the flows in 
semidetached binaries is examined in only a few of 
these, since many of these studies were performed 
in a restricted formulation in which the process of 
establishment of flow formation was not taken into 
consideration. The variety of binary systems studied 
and approaches used makes it impossible to general- 
ize the results obtained, or to make at least qualita- 
tive predictions about the flow structures expected for 
semidetached systems. At the same time, such studies 
are extremely important, since data on the structure 
of gas flows are necessary for interpretation of the ob- 
servational data, and the time-consuming nature of 
3D modeling makes it impossible to compute the flow 
pattern for each particular system. 

Here, we consider the results of 3D numerical mod- 
cling of the flow of an ideal gas in a semidetached bi- 
nary system assuming the absence of magnetic field 
and radiative heating and cooling of the gas. To make 
the study more systematic, we generalize the results 
of eight of our own calculations, and have also used, in 
some cases, the results of other studies. We identified 



the basic parameters of this problem (Section 2) and 
successively examined the effect of each parameter on 
the flow structure (Section 3). Features of the flow 
brought about by the influence of various parameters 
and the typical characteristics of the flow structure 
are described in Section 4. 

2. Parameters of the Problem 

In order to perform gas-dynamical studies of the 
flow structures in binary systems, it is necessary to 
specify values of the parameters describing the system 
under consideration and of the parameters used in the 
mathematical model used. 

2.1. Parameters of binary system 

Numerical analyses of the matter flows in binary 
systems to not require the specification of all seven 
orbital elements for the binary star, since the orbital 
plane is assumed to be known in advance. In this 
case, we can describe the motion of the system's com- 
ponents using the mass of the donor star Mi and the 
mass of the accretor star Mi (or equivalently, the to- 
tal mass of the system M and the mass ratio of the 
components q — M2/M1), the semimajor axis of the 
orbit A, the rotation period P or i, (or, equivalently, 
the angular velocity £1 = 27r/P or b), and the orbital 
eccentricity e. Note that the application of Kepler's 
third law decreases the required number of param- 
eters, making it possible to use, instead of the four 
parameters M, q, A, and Cl, any three of these to 
describe the system in a point-mass approximation. 
The geometric sizes of the components should also be 
included among the parameters describing the sys- 
tem. For the semidetached systems considered here, 
we assume that the donor star fills its Roche lobe and, 
therefore, its geometric characteristics are completely 
specified by the parameters of the system. The size 
of the accretor star i? 2 is an independent parameter 
that must be given separately. 

2.2. Model parameters 

In our numerical model, we use the Euler equa- 
tions in the form of integral conservation laws to de- 
scribe the flow of a nonviscous, non-heat-conducting 
gas in a binary system without magnetic or radiation 
fields. This permits us to obtain discontinuous (non- 
isentropic) solutions satisfying the Rankine-Hugoniot 
conditions at the discontinuities (shock waves and 
contact discontinuities). We complete the system of 
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gas-dynamical equations with the equation of state 
of an ideal gas with adiabatic index 7: P = (7 — l)pe, 
where P is the pressure, p is the density, and e is the 
specific internal energy. 

To solve the system of gas-dynamical equations, 
we must specify the boundary conditions at the sur- 
faces of the system's components and at the surface of 
the computation region. The states of the gas in the 
near-surface layers of the donor and accretor stars are 
determined by the gas densities p\ and p 2 , the temper- 
atures Ti and T 2 (or, equivalently, the sound speeds 
ci and c 2 ), and the gas- velocity vectors vi and v 2 . 
The boundary conditions at the stellar surfaces are 
obtained by the standard procedure of solving Rie- 
mann problem between the states immediately un- 
der the stellar surface and the first computation node 
(see, for example, [2, 3]). A free-outflow condition is 
specified at the outer boundary of the computation 
region. 

Thus, in our formulation, the motion of a nonvis- 
cous gas in a binary system without magnetic and 
radiation fields depends on five parameters specifying 
the binary system: q, A, CI, e, R2 and seven model 
parameters: 7, pi, c\, vi, pi, c 2 , v 2 . 

Understanding the necessity of studying succes- 
sively the roles of individual parameters, we introduce 
a number of additional constraints that decrease the 
total number of parameters. In particular, we assume 
that we can neglect the stellar wind produced by the 
accretor star. This is standard in studies of flows in 
semidetached binaries, and is widely used in computa- 
tions (see, for example, [12-19]). In this case, the only 
influence of the accretor on the flow of matter is grav- 
itational, and we can adopt a condition of complete 
accretion at its surface, i.e., the condition that all gas 
reaching the surface of this component is absorbed. 
In order to be able to perform a stage-by-stage inves- 
tigation of the effect of various parameters on the flow 
structure, we also add a number of constraints on the 
type of system under consideration: we assume that 
(i) the orbits of the components are circular, which al- 
lows us to fix the parameter e = 0; (ii) the accretor is a 
compact object (a white dwarf, neutron star, or black 
hole), so that, when i? 2 *C A, the general structure of 
the flow is determined only by the gravitational field 
of the accretor and, accordingly, the parameter R 2 
can be excluded; (iii) the rotation of the donor star is 
synchronized with the orbital motion of the system, 
so that, in a reference frame rotating with the angular 
velocity of the system, the gas velocity at the surface 



of the donor star is directed along the normal to this 
surface, and we can replace the vector parameter vi 
by the scalar parameter v\ (equal to the velocity com- 
ponent perpendicular to the donor-star surface). 

In the above formulation, the flow structure in the 
binary system depends on seven parameters: q, A, 
Ct, 7, pi, ci, v\. It is noteworthy that the solution 
is not affected by the boundary density at the donor 
star, since the system of equations can be scaled with 
respect to p (and simultaneously with respect to the 
pressure P). In the numerical model, an arbitrary 
rho\ value can be used. When considering a specific 
system with known mass-loss rate, however, we must 
scale the calculated density values by the ratio of the 
true and model densities at the surface of the donor 
star if we wish to find the real density values for the 
system. For convenience in the numerical solution, 
the original system of equations can be written in di- 
mensionless form. Following the usual practice (see, 
for example, [2, 3]), we use the distance A between 
the components and the reciprocal of the angular ro- 
tation velocity of the system O -1 , to make the equa- 
tions dimensionless with respect to space and time, 
respectively. In this case, the remaining six parame- 
ters determining the flow structure can be reduced 
to a combination of four dimensionless parameters 
q = M 2 /M 1 , e = d/(ClA), M = Vi/ci, 7, and so- 
lutions obtained for two systems with the same pa- 
rameter values will differ only by a scaling factor. 

3. Results 

To determine the influence of various parameters 
on the flow structure, we carried out 3D numerical 
modeling of several binary systems. We successively 
varied each of the parameters while fixing the values 
of the remaining parameters. It is evident that in 
a full analysis, the number of values considered for 
each parameter (in the total characteristic range of 
values) must be as large as possible. Unfortunately, 
3D modeling requires substantial computational re- 
sources, and we restricted ourselves to only eight cal- 
culations. The parameters used in these calculations 
are summarized in Table I. The limited number of 
calculations makes it impossible to study the depen- 
dence of the results on a selected parameter in detail. 
Nevertheless, analysis of our results enables identifica- 
tion of the parameters most strongly affecting the flow 
structure, and determination at a qualitative level of 
trends in the solution behavior. 
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X/A 

Fig. 1. — (a) density contours and velocity vectors in the equatorial plane of the binary system for calculation 
1 specified in Table I. The vector in the upper right corner corresponds to the velocity AA ■ f2. Three Bowlines, 
denoted 'a', and 'c' illustrate the directions of the flow of matter in the system. The dashed lines are the Roche 
equipotential surfaces. 



Taking into account the computer-intensive nature 
of the problem, we included our previous results [22- 
25] in our analysis. Additional computations were 
made in three cases (computations 3, 4, 5 in Table I), 
in order to fill out the field of parameter values. We 
also used results presented by other authors when 
their parameters differed from those in Table I. 

The computation results for each of the eight cases 
presented in Table 1 are given in Figures 1-8. The 
upper panels (Figures la-8a) depict density contours 
and velocity vectors in the equatorial plane of the 
system illustrating the morphology of gas flows in the 
system. A vector showing the velocity scale is given 
in the upper right corner of each plot. In all cases, the 
length of this vector corresponds to the velocity AA-Q. 
We also show in the figures three flowlines marked by 
letters 'a', l b\ and 'c', which illustrate the directions 
of the flow of matter in the system. Density contours 
and velocity vectors in the plane passing through the 



Table I. 



Case N 


q 


e 


M 


7 


1 


l 


0.15 


0.083 


1.2 


2 


l 


0.15 


0.083 


1.01 


3 


l 


0.016 


0.083 


1.2 


4 


l 


0.016 


0.083 


1.01 


5 


l 


0.016 


1 


1.01 


6 


5 


0.016 





1.01 


7 


5 


0.016 


1 


1.01 


8 


5 


0.04 


1 


1.01 



accretor and perpendicular to the line connecting the 
centers of the two components are shown in the lower 
panels (Figures lb-8b). The velocity scale is specified 
by the vector in the upper right corner of these panels, 
which corresponds to the velocity 2A-VL. For conve- 
nience, the range of density variation in Figures lb- 
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Fig. 1 - continued. — (b) density contours and velocity vectors in the plane passing through the accretor star and 
perpendicular to the line connecting the centers of the system's components. The vector in the upper right corner 
corresponds to the velocity 2A ■ il. The range of the density variations is 0.1 - 0.0003/?o- The thick solid line 
corresponds to the density O.OOlpo- 



8b is chosen to be the same, and to be equal to 0.1 
- 0.0003po- The thick line corresponds to the density 
O.OOlpo- The flow pattern in the figures is stationary 
and corresponds to a steady-state flow regime. 

Analysis of the results presented in this paper en- 
ables investigation of the influence of various param- 
eters on the flow structures in semidetached binary 
systems. 

3.1. The effect of 7 

The dependence of the flow pattern on the adopted 
value of the parameter 7 can be established by com- 
paring the results of calculations in cases 1 and 2, and 
in cases 3 and 4. We varied only the parameter 7 in 
these pairs of cases, with all other parameters having 
fixed values. 

Analysis of the results presented in the figures 
shows that, in the quasi-stationary case (0.001 po), 
the flow pattern is characterized by the presence of 
an accretion disk and circumbinary envelope (Fig- 



ures 2 and 4). The details of the flow structure in 
the cases under consideration are largely determined 
by the presence of the circumbinary envelope, and 
were considered by us in [22-25]. Below, we give 
only a brief description of the main features of the 
flow for this case, referring the reader to [22-25] for 
a more detailed discussion. Analysis of the gas flows 
allows us to divide the matter in the stream flowing 
from L\ into three parts: the first produces a quasi- 
elliptical accretion disk (flowline 'a') and continues 
to participate in the accretion process, losing angular 
momentum under the action of viscosity; the second 
(flowline '&') flows around the accretor outside the 
disk; the third part of the stream moves toward the 
Lagrange point L2 (flowline 'c'), and a significant por- 
tion of this matter changes the direction of its motion 
under the action of the Coriolis force and remains in 
the system. The matter remaining in the system and 
not involved directly in the accretion process (flow- 
lines '6' and 'c') forms the circumbinary envelope of 
the system. A significant portion of the envelope gas 



5 




0.3 
0.2 
0.1 

< 

^> 0.0 



-0.4 



\ \ \ \ \ i r t r / / / 
\\\ \ \ i i t r / // // 



\ \ \ \ \ \ 

v v \ \ \ A 1 

^ v v \ \ \ 1 

. ^ v s \ T 




r / ///////•/ 
/ / / / / / / / S / 
r / r / / s / s y / 

7 / / / S S / S S s 

r / s / S s» 

s* ^ — * ^> — » 



— -. — * — -» 
\ \ S N. S ^» ^> "--» 

S / / / / 1 I i l \ N 



•s s s / / / / J 
s s / J I I J 



I 

-0.2 



0.0 

Y/A 



0.2 



0.4 



Fig. 3. — Same as Figure 1 for calculation 3 from Table I. 
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Fig. 4. — Same as Figure 1 for calculation 4 from Table I. 
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Fig. 5. — Same as Figure 1 for calculation 5 from Table I. 
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Fig. 6. — Same as Figure 1 for calculation 6 from Table I. 
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Fig. 7. — Same as Figure 1 for calculation 7 from Table I. 
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Fig. 8. — Same as Figure 1 for calculation 8 from Table I. 
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interacts with the matter flowing from the surface of 
the donor star, and substantially changes the mass ex- 
change regime in the system [23, 25]. Another part of 
the circumbinary envelope (flowline '6') passes around 
the accretor and undergoes a shock interaction with 
the edge of the jet facing the direction of the orbital 
motion. This interaction leads to the absence of a 
shock interaction of the stream flowing from L\ with 
the gas of the forming disk. The stream, deflected 
by the envelope gas, comes at the disk along a tan- 
gent, and does not produce a shock disturbance at 
the disk edge ("hot spot"). At the same time, the 
interaction of the stream with the circumbinary en- 
velope results in the formation of an extended shock 
wave (denoted I in the figures) along the edge of the 
stream, which can be considered equivalent to a " hot 
spot" on the disk in terms of its observational mani- 
festation [22, 24]. The interaction of the circumbinary 
envelope with the material in the stream also leads to 
the formation of shock wave II. 

The flow pattern undergoes qualitative changes in 
calculations with the higher index 7 = 1.2 (Figures 1 
and 3). This is revealed primarily in the absence of 
the accretion disk. Analysis of the flowlines presented 
in Figures la and 3a shows that some portion of the 
stream material (flowline 'a') flows directly onto the 
accretor, whereas another portion (flowline '&') passes 
around the accretor and undergoes a shock interaction 
with the stream flowing from Li, without forming a 
disk. This fact is supported by the appearance of the 
density and velocity-vector fields in the YZ plane in 
Figures lb and 3b. Comparison of these with Fig- 
ures 2b and 4b provides conclusive evidence for the 
absence of an accretion disk for the higher 7 value. 
Note that this result is in agreement with [14, 16], 
which demonstrated that a disk forms only when 7 
is close to unity (7 < 1.1). Additional distinctions in 
the flow pattern for the different 7 reveal themselves 
in changes of the directions of flows in the system. In 
particular, analysis of the gas motions in the Z di- 
rection shows that, for the higher 7, the gas moves 
predominantly away from the equatorial plane (Fig- 
ures lb and 3b), whereas for 7 ~ 1, it moves from 
higher heights toward the accretion disk (Figures 2b 
and 4b). 

Along with the qualitative differences in the flow 
pattern, the variations in the value of 7 give rise to 
significant quantitative changes. In particular, for the 
low value of 7, the influence of the circumbinary en- 
velope on the surface of the donor star substantially 



increases. This leads to an increase in the rate of mass 
loss by the donor star and, as a consequence, to an 
increase in the matter accretion rate by a factor of 
1.5-2. 

3.2. The effect of M 

Pairs of calculations that can be used to assess the 
role of M. are denoted in Table I by the numbers 4, 
5 and 6, 7. To study the role of the parameter A4, 
we compared calculations in which the gas velocity at 
the near-surface layer of the donor star was taken to 
be small (the Mach number M. = 0.08 in calculation 
4 and the Mach number M = in calculation 6) 
with calculations in which the gas velocity component 
normal to the surface was equal to the local sound 
speed (A4 = 1 for models 5 and 7 in Table I). 

The results indicate that a number of minor quanti- 
tative changes arise in the flow pattern. In particular, 
as the Mach number increases, we observe a growth 
in the rate of mass loss by the donor star due to the 
" stripping" of its stellar atmosphere by the envelope 
gas. This effect is quite natural, since an increase in 
the near-surface gas velocity is accompanied by an in- 
crease in the overall extent of the atmosphere, which 
leads to more intense stripping of atmospheric mate- 
rial by the gas of the circumbinary envelope. Note, 
however, that the growth of the mass-loss rate of the 
donor star does not influence the accretion rate, which 
remains virtually constant as the Mach number varies. 

At the same time, the general flow pattern retains 
all the previously determined qualitative character- 
istics (see Figures 4 and 5 and also Figures 6 and 
7). Moreover, the basic parameters of the gas flows, 
dimensions and position of the accretion disk, and po- 
sition and intensity of shock waves / and II remain 
virtually unchanged. This, in turn, suggests that M 
only weakly influences the flow structure. 

3.3. The effect of e 

The parameter e = c\/(ClA) was first introduced 
by Lubow and Shu [26], and was used as a small 
parameter. Let us estimate the characteristic range 
of this parameter for various types of binary sys- 
tems. Combining the expressions QA = yj GM/A 
and ci = y/-fRTi , we can represent e as 

e = O.OOST^P^M^V 72 • 

Here, the orbital period P is expressed in years and 
the total mass M of the system in solar masses M Q . 
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Fig. 9. — One-dimensional density profiles for cal- 
culations 2 and 4 specified in Table I. The density 
values are taken along the straight line in the equa- 
torial plane and passing through the point 0.6A per- 
pendicular to the line connecting the centers of the 
components of the system. 

Let us assume that the gas temperature at the donor 
star is 10 4 K. The range of e values characteristic of 
this case is presented in Table II for various P and 
M. Analysis of the data in Table II shows that the 
range of e is rather small (~ 0.001 -j- 0.3). However, e 
cannot always be considered a small parameter. 
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M = 


Mq 


M = 1OM 


M = 1OOM 
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year 


0.3 


year 


0.13 


year 


0.065 


month 


0.13 


month 


0.05 


month 


0.028 


day 


0.04 


day 


0.02 


day 


0.009 


hour 


0.014 


hour 


0.006 


hour 


0.003 


minute 


0.0037 


minute 


0.0017 


minute 


0.0008 



To determine the effect of the parameter e, we com- 
pare calculations 2 and 4, and 7 and 8. In the first 
pair of calculations, we varied e from 0.15 to 0.016, 
and in the second pair, from 0.016 to 0.04. Analysis 
of the results (Figures 2, 4, 7, and 8) shows that the 
qualitative flow pattern remains unchanged for vari- 
ous e: an accretion disk forms in the system; there 
is no "hot spot" in the region of interaction of the 
stream with the disk; the directions of the flows coin- 
cide; and the forming circumbinary envelope interacts 



with the stream, producing the shock waves / and II. 
In addition, a number of quantitative changes are ob- 
served as e is increased. In particular, the cross sec- 
tion of the stream increases in the calculations with 
the larger e. This is supported by Figure 9, which 
presents one-dimensional density profiles for calcula- 
tions 2 and 4. The density values are taken along the 
straight line in the equatorial plane passing through 
the point 0.6A perpendicular to the line connecting 
the centers of the system's components. This conclu- 
sion is consistent with the results of Lubow and Shu 
[26], who showed that, for small e, the width of the 
stream is e. In addition to the increase in the stream 
size for the larger value of e, we note an increase of 
the disk thickness and the accretion rate. 

3.4. The effect of q 

In this study, we examined systems in which the 
parameter q — M 2 /Mi takes on only two values: 
q = 1 in calculations 1-5 and q = 5 in calculations 
6-8. Assessment of the role of this parameter can 
be made using calculations 5 and 7, in which all pa- 
rameters except q were the same. Comparison of the 
results presented in Figures 5 and 7 shows that the 
qualitative flow pattern remains unchanged, although 
some quantitative variations are present. This indi- 
cates that the calculation results depend only slightly 
on the adopted ratio of the masses of the components 
when q > 1. 

The dependence of the flow pattern on q when 
q < 1 requires additional study. According to [27, 
28], when q < 1, the flow pattern undergoes signifi- 
cant changes, which is reflected, in particular, in clear 
observational manifestations. In addition, consider- 
able tidal disturbances are possible when the mass of 
the donor star is large. A detailed study of the influ- 
ence of the parameter q will be presented in a future 
paper. 

4. Conclusion 

Analysis of the 3D gas-dynamical calculations per- 
formed here enables us to sort the parameters of the 
problem according to their influence on the flow pat- 
tern in a semidetached binary system. The variation 
in the adiabatic index 7 has the largest effect on the 
solution. The influence of M and e on the flow pat- 
tern is small, and leads only to unimportant quan- 
titative variations of the solutions without changing 
the overall qualitative flow patter. The ratio q of the 
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masses of the two components also does not exert a 
significant influence for the region q > 1. 

Since qualitative changes of the flow pattern de- 
pend only on the adopted value of the parameter 7, 
let us consider the mechanism of its action upon the 
flow and summarize the main features of the gas-flow 
structure that should be observed in semidetached bi- 
naries. We considered the flow of a nonviscous, non- 
heat-conducting gas in a binary system without mag- 
netic and radiation fields. In the model, the assump- 
tion that the gas is ideal formally implies neglect of 
radiative losses, although small values of the adiabatic 
index 7 imitate energy losses in the system (see, for 
example, [3]). In particular, when 7 ~ 1, a solution 
close to the isothermic solution is realized [28]. In 
the above calculations, we used the two 7 values 1.2 
and 1.01, and, thus, obtained solutions corresponding 
to models with different energy-loss rates. The flow 
patterns calculated in these two cases display funda- 
mental differences. For example, in the calculations 
performed with 7 = 1.2, the accretion disk does not 
form at all. The increase in energy losses in the sys- 
tem in the calculations for the close to isothermal case 
(7 <~ 1) leads to the formation of an accretion disk and 
a corresponding change in the accretion mechanism. 

In typical binaries, the gas between the compo- 
nents is ionized, and the corresponding energy-loss 
mechanisms work efficiently, so the real flow pattern 
should be described by a close to isothermal model 
[30, 31]. As we noted above, in this case, the flow 
pattern is characterized by the following features: the 
stream of gas flowing from L\ forms both an accretion 
disk and a circumbinary envelope in the system. Fur- 
ther, the envelope gas exerts a considerable influence 
on the structure of gas flows in the system [23-25] . In 
particular, the envelope gas interacts with the stream 
of matter flowing from the vicinity of L\ and deviates 
it, giving rise to a shockless (tangential) interaction 
of the stream with the outer edge of the forming ac- 
cretion disk and, therefore, to the absence of a "hot 
spot" in the disk. At the same time, the interaction 
of the envelope gas with the stream is responsible for 
the formation of an extended shock wave at the edge 
of the stream, whose observational manifestations can 
be considered equivalent to those of a hot spot in the 
disk [24]. 
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